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It.  AstraAcr 


The  characterization  of  ionospheric  scintillations  in  terms  of  scintillation 
Index  is  discussed.  Scintillation  Index  cannot  be  directly  applied  to  the  design 
of  systems  that  use  transionospheric  propagation.  The  relationships  of  scintil¬ 
lation  index  to  system  design  parameters,  cumulative  amplitude  distributions, 
and  fade  levels,  are  shown.  The  determination  of  a  spectral  index  is  explained 
and  its  use  for  extrapolating  the  results  of  ionospheric  scintillation  measure¬ 
ments  to  a  desired  frequency  Is  shown. 
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Notes  on  the  Relationship  of  Scintillation  Index  to 
Probability  Distributions  and  Their  Uses  for  System  Design 


I.  INTRODUCTION 

Ionospheric  sclntllUtions  can  be  a  serious  problem  to  many  systems  that 
transmit  through  the  ionosphere  and  are  sensitive  to  fading  of  the  signal.  Some 
communication  systems  are  designed  with  as  little  as  a  4-dB  margin  and  scintilla' 
tlons  can  often  cause  this  level  to  be  esceeded  even  at  mid'latitudes  for  VHP  and 
UHF  transmissions.  Equatorial  and  auroral  stations  have  shown  greater  changes 
at  even  higher  frequencies.  Much  data  that  has  been  taken  for  characterizing 
scintillations  is  not  in  a  form  which  can  be  directly  used  by  the  engineer  for 
system  design  or  the  physicist  interested  in  mechanisms  governing  the  generation 
of  scintillations.  It  is  th  '  intent  of  this  note  to  show  the  relationships  between 
scintillation  indices  and  amplitude  distributions  and  ways  in  which  existing  data 
can  be  made  useful  for  systems'  application.  Since  a  given  application  will  often 
be  for  a  frequency  which  is  different  than  that  described  by  the  data,  a  method  is 
presented  for  scaling  the  data  to  other  frequencies. 


(Received  for  publication  2  January  1974) 
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2.  SCINTILLATION  INDEX 


The  morphology  of  Ionospheric  sclntlUiitlons  has  been  usefully  described  by 

scintillation  indices  that  are  a  measure  of  the  fluctuations  Imposed  on  a  signal  as  It 

traverses  the  Ionosphere.  A  simple  index  was  adopted  by  AFCRL  for  uniformly 

1  2 

scaling  chart  recordings  of  the  strength  of  satellite  beacons.  ‘  Scintillation  index 

(Sn  was  defined  as  where  is  the  power 

,  max  min  max  min  max 

amplitude  of  the  third  peak  down  from  the  maximum  excursion  of  the  scintillations 

and  is  the  power  amplitude  of  the  third  level  up  from  the  minimum  excursion 

in  the  sample  period,  typically  15  minutes.  It  was  found  that  scaling  the  chart 

records  was  considerably  simplified  by  just  measuring  the  decibel  change  between 

P  and  P  ,  .  By  assuming  equal  percentage  changes  from  the  average  level  for 
max  min 

^max  ^min'  *  relationship  between  dB  and  percentage  canbe  determined.  This 
is  shown  in  Table  1  and  Figure  1. 


Table  1.  Conversion  Table;  _  -  P,^,„(dB)  vs  Si  (%) 

max  min 


-dB 

♦  dB 

P  -  P  ,  (dB) 

max  min 

SI  {%) 

.09 

.09 

.  18 

2 

.  22 

.21 

.43 

5 

.46 

.41 

.87 

10 

.71 

.61 

1.32 

15 

.97 

.79 

1.76 

20 

1.29 

.97 

2.22 

25 

1.  55 

1.  14 

2.69 

30 

1.87 

1.30 

3.17 

35 

2.22 

1.46 

3.68 

40 

2.60 

1.61 

4.21 

45 

3.01 

1.76 

4.77 

50 

3.47 

1.90 

5.37 

55 

3.98 

2.01 

5.99 

60 

4.56 

2.  18 

6.74 

65 

5.23 

2.?1 

7.54 

70 

6.02 

2.43 

8.45 

75 

6.99 

2.55 

9.54 

80 

8.24 

2.87 

10.91 

85 

10.00 

2.79 

12.79 

90 

13.00 

2.  90 

15.  90 

95 

20.00 

2.99 

22.99 

99 

30.00 

3.00 

33.00 

99.9 

1.  Whitney.  H.  E.  and  Mslik,  C.  (1868)  A  Proposed  Index  for  Measuring 

Ionospheric  Scintillation.  AFCRL*68-0138. 

2.  Whitney.  H,E..  Aarons,  J.  and  Malik,  C.  (1969)  A  proposed  Index  for 

measuring  ionospheric  scintillations.  Planet.  Space  Scl. .  17:1069-1073, 


Figure  I,  Graph  for  the  Converaion  of  Pmax  *  ^mln  Scintillation  Index. 
PjTiax  *  **min  P®®^  P*®**  excursion  of  a  scintillating  signal  and  is 

measured  In  decibels  based  on  an  amplitude  calibration  of  the  chart  record 


3.  CUMULATIVi;  AMPLITUDE  PROBABILITY  DISTRIBUTION  FUNCTION 

While  the  characterization  of  ionospheric  scintillations  has  been  extremely 
useful  for  describing  the  morphology  of  the  Irregularity  structure,  it  has  not  des¬ 
cribed  the  fading  characteristics  In  a  form  suitable  for  engineering  design,  such 
as  the  fading  margin  necessary  to  overcome  scintillations  as  a  function  of  local 
time,  latitude,  season,  magnetic  activity,  and  so  forth.  The  cumulative  amplitude 
probability  distribution  function  (cdf)  Is  a  form  of  data  that  Is  readily  applied  to 
system  design.  The  relationship  between  SI  and  cdf  has  been  determined  for  a 
large  number  of  experimental  distributions  of  synchronous  satellite  signals  at 
136  MHz.®  General  agreement  between  the  experimental  distributions  and  the 
Nakagami  m  distribution  was  found;  examples  are  shown  In  Figure  2.  Figure  2 

3.  Whitney,  H.E.,  Aarons,  J. ,  Allen,  R.S.  and  Seem ann,  D.R.  (1972)  Estima¬ 
tion  of  the  cumulative  amplitude  probability  distribution  function  of  Iono¬ 
spheric  scintillations.  Radio  Science,  7:1095-1104. 
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also  lUustr.tes  the  scaling  of  scintillations  for  values.  The  punched 

card  format  for  computer  processing  Is  given  In  Appendix  A.  Figure  3  shows  the 
theoret'cal  curves  for  the  Nakagaml  distribution  over  a  range  of  m  values.  If  ex- 
perlme  tb'  data  does  not  cover  sufficient  amplitude  range.  Figure  3  can  be  used  to 
estimate  xtreme  values.  Figure  4  shows  the  grouping  of  scintillation  index  values 


9999  99  90  50  10  5  I  0.1  OOl 


Figure  3.  Theoretical  Nakagaml  m  Distribution 
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poKenr  Of  time  sicnal  amalituoe  exceeoco  ohoinate 

Figure  4.  Model  Distributions  for  Ionospheric  Scintillations 
at  136  MHz 

that  was  used  to  facilitate  statistical  studies  and  to  generate  the  model  distributions. 
The  curves  are  median  distributions  for  each  SI  group.  The  solid  corves  denote 
the  range  of  the  experimental  data  and  the  dashed  curves  denote  an  extrapolation 
based  on  the  Nakagami  distribution  for  the  value  of  m  in  parenthesis.  The  standard 

3 

error  of  the  mean  equals  (7  Additional  details  are  given  in  Whitney  et  al. 

Our  analysis  of  experimentally  obtained  amplitude  distributions  at  137  MHz 
has  shown  that  there  is  close  agreement  with  the  Nakagami  m  distribution.  It  has 
also  been  shown  that  =  l/y  m  for  m  s  0.  5,  where  is  the  notation  used  by 


4,  Chytil,  B.  (1967)  The  distribution  of  amplitude  scintillation  and  the  conversion 
of  scintillation  indices.  J,  Atmos,  Terr.  Phys. ,  29:1175-1177. 
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Figure  5.  The  Relation  Between  and  m 


4.  RELATIONSHIP  OF  Si  TO  m 

Fifty -two  15-mln  samples  of  scintillation  data  were  processed  by  the  com¬ 
puter  to  obtain  a  cdf  and  the  best  fit  to  a  particular  m  value.  The  samples  were 
also  scaled  for  SI  (dB).  A  smoothed  curve  giving  the  relationship  of  scintillation 
Index  (dB)  to  the  Nakagaml  m  parameter  is  shown  in  Figure  6. 

Since  Figure  1  shows  the  relationship  between  SI(dB)  and  SI(7i)>.  Figure  6  can 
be  replotted  to  give  the  variation  of  SI(7«)  with  m  as  shown  in  Figure  7. 


S.  Briggs,  B.H.  and  Parkin,  I.  A.  (1963)  On  the  variation  of  radio  itar  and  satel¬ 
lite  scintillations  with  zenith  angle,  J.  Atmos,  Terr,  Phys. ,  35:339-366. 
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1  \  2  4««I0  20  404OMI0O  200  4 

m 

Figure  6.  Variation  of  Nakagami  m  Parameter  with  Scintillation  Index,  m  was 
determined  by  computer  processing  of  137 -MHz  data;  SI  was  obtained  by  sealing 
strip  charts  for  fifty-two  15 -min  samples 


I  2  4  «  10  20  40  ro  100  200  SOO 

m 

Figure  7.  Scintillation  Index  vs  m  Taken  from  Figure  1  and  Figure  6 


S.  RELATIONSHIP  OF  SI  TO  PERCENTILE  POINTS 


Figure  8  shows  a  histogram  of  the  percentile  points  of  the  fifty-two  amplitude 
distributions  that  had  a  dB  change  equal  to  P  -  P_,  .  Each  of  the  52  samples 
used  for  Figure  6  had  an  St  value  and  an  amplitude  distribution.  The  percentile 
points  (plus  and  minus  from  the  median  level)  were  recorded  from  the  distribution 
that  gave  a  dB  change  equal  to  the  scaled  value  from  the  chart  records.  It  shows 
that.  In  general,  the  value  of  SI  using  the  third  peak  definition  agrees  with  the  one 
percent  points  on  the  cdf. 


Figure  8.  Histogram  of  the  Amplitude  Distribution  Percentile 
Points  that  Correspond  to  the  Pmax  ~  Pmin  ^^B)  Levels  for 
the  Fifty-Two  15-mln  Samples  in  Figure  6 


6.  RELATIONSHIP  OF  .SI  TO  FADE  LEVELS 


The  conversion  between  Sl(7o)  and  SI<dB),  that  is  P_  -  P_.  ,  was  based  on 

max  min 

equal  percentage  changes  from  the  average  level  (Table  1).  Since  the  negative 
changes  or  fade  levels  were  recorded  for  the  data  for  Figure  8,  we  can  check  the 
experimentally  determined  change  against  the  assumption  of  equal  percentage  change. 
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Th«  results  ere  shown  in  Figure  9  where  the  calculated  fade  level  Is  compared  with 
the  measured  fade  level  for  various  SHdBl  levels.  It  Illustrates  that  anSI  >  8.  5  dB 
(75  percent)  would  give  a  6-dB  fade  according  to  the  calculated  relationship  (solid 
line),  but  the  experimental  data  would  require  an  SI  >  10  dB  (82  percent)  to  give  a 
S-dB  fade. 


Figure  9.  Comparison  of  the  Calculated*  with  the  Measured-  Fading 
Range  as  a  Function  of  Scintillation  Index 


7.  FREQUENCY  DEPENDENCE 

It  is  important  to  be  able  to  extrapolate  the  results  of  scintillation  studies  at  a 
specific  frequency  to  another  frequency  of  interest.  In  order  to  accomplish  this, 
the  frequency  dependence  or  spectral  index  of  scintillations  must  be  known.  Bared 

on  scintillation  index  S,  a  spectral  index  can  be  expressed  as  n  >  (log  (S, /S.,)]/(log 

1/2  ^  * 

Since  S  >  (1/m)  ,  a  spectral  index  in  terms  of  the  m  parameter  can  be 

obtained  ast}^  <  [log  (m  j/m^ll/llog  (fj/f^)].  These  expressions  are  valid  as  Long 

as  scintillation  index  and  m  are  a  constant  power -law  function  of  frequency. 


14 


Satellites  ATS-3  and  ATS-S  had  transmissions  at  137  and  412  MHt  which 
allowed  the  measurements  of  an  m  value  for  lid  simultaneous  15>mln  amplitude 
distributions  at  each  frequency.  The  spectral  index  was  then  calculated  and  the 
results  are  shown  in  the  histogram.  Figure  10. 


Figure  10.  Histogram  of  Spectral  Index 


a  USES 

ai  Extrapolation  of  the  EffacU  of  SdntilUtioM  to  a  De^ed  Frequency 

Suppose,  fcr  example,  a  short  sample  <15  min)  of  scintillation  data  is  measured 
at  137  MHx  to  have  a  cdf  with  an  m  «  1. 5  and  it  desired  to  estimate  the  cdf  at  254 
MHz.  Under  the  assumption  that  the  measured  spectral  index  ^  is  a  constant 
power  law  function  of  frequency,  we  can  use  the  average  value  of  2. 62  from  Figure 
10  to  evaluate  the  following  expression; 


IS 


in„ 

m. 


'  m 


log 


‘H 


^  ?“ 


IHu 

log  ^  c  2.62(.268)  -  .703 

"^L 


rtijj  -  7.  57 


■■  254  MH2 
*  137  MHz 

Im  =  2.62 


The  resulting  distribution  for  the  calculated  m  can  then  be  read  or  interpolated 
from  the  theoretical  m  curves  shown  in  Figure  3. 


8.2  Convenioo  of  Occurrenee  Statirtiei  of  SI  (o  a  edf 

We  have  measured  137  MHz  scintillations  at  Narssarssuaq.  Greenland,  for 
several  years  and  have  selected  the  following  statistics  to  show  the  technique  for 

3 

converting  to  a  cdf.  The  technique  has  been  shown  to  give  sufficient  accuracy  for 
engineering  purposes  by  the  use  of  the  models  which  are  determined  for  Hamilton 
data  and  shown  in  Figure  4. 

Table  2.  Occurrence  for  Narssarssuaq;  February-April  1968*1972 


K  =  0  -  3 

2200 

-  0200  LT 

SI  Group 

SI% 

%  Occurrence 

0 

5.31 

1 

8.85 

2 

12.48 

3 

60-79 

16.51 

4 

80-86 

11.66 

5 

^90 

45.  18 

The  percent  occurrence  of  the  15-mln  scintillation  Indices  are  given  for  the  4-hour 
local  time  block  and  a  two-month  period  for  the  Indicated  years  under  quiet  mag¬ 
netic  conditions.  To  determine  a  cdf  which  represents  this  sorting  of  SI,  the 
percent  occurrence  of  each  SI  group  Is  multiplied  by  its  probability  at  each  dB  level 


16 


M  given  by  the  models  In  Figure  4.  end  the  values  ere  summed  to  give  s  composite 
curve.  Expressed  methemeticelly: 


cdf(dB) 


Z 


Croups  0 


r(G)P(dB). 


where  f(C)  is  the  frequency  of  occurrence  ofe  group  end  P(dB)  is  the  probebility  or 
percent  at  time  et  each  dB  level  from  Figure  4.  The  dashed  curve  in  Figure  1 1 
then  gives  the  distribution  for  the  occurrence  of  SI  groups  listed  above. 


PERCENT  OF  TIME  SIONAL  AMPLITUDE  EX9EEOED  ORDINATE 

Figure  11.  Extrapolated  Scintillation  Amplitude  Distribution 
Based  on  Measured  Data  from  Narssarssuaq  at  137  MHz 


8.3  Chw  dM  Omaiiiaw  SutiMici  for  137-MHi  SI  ••  LMod  bi  PieetdiBK  Soctioa,  FM 
Expoctad  DNti'fcuMue  at  SS4  MHs  for  the  StatMea 

Again,  the  assumption  is  made  that  the  spectral  index  is  a  constant  power  law 
function  of  frequency  and  is  given  by  the  average  value  from  Figure  10,  that  is, 

1}  s  2.62.  In  the  same  manner  that  an  individual  distribution  was  converted  to  a 
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different  frequency  es  In  Section  8. 1.  then  ell  the  model  distributions  of  Figure  4 
can  be  converted  to  their  equivalents  at  254  MHi  by  *  5. 05  ^  8  The 

m  values  for  each  SI  group  are  listed  in  Table  3  for  137  MH-  id  the  converted 
value  to  254  MHz. 


Table  3.  The  m  Values  for  Each  SI  Group 


Group 

“"l 

(13?  MHz) 

*"h 

(254  MHz) 

1 

40 

200 

(negligible) 

2 

IS 

75 

3 

10 

SO 

4 

4 

20 

5 

2 

10 

The  percent  of  time  values  for  each  dB  level  for  each  curve  were  read  and 
then  multiplied  by  the  percent  occurrence  of  the  S!  groups  to  give  the  resulting 
distribution  at  254  MHz  shown  by  the  solid  curve  in  Figure  11. 

8.4  Given  an  Sa  at  400  MHz  (Lineola  Laboratery  Data),  Find  the  Equivalent 
AFCRLSIal254MHt 

2 

Given  an  S^  at  400  MHz,  we  can  calculate  an  m  at  400  MHz  from  m  *  l/S^  or 
Figure  S.  Assuming  that  the  frequency  dependence  experimentally  determined 
from  137 -MHz  and  4 12 -MHz  data  is  valid,  the  relationship  of  m^.,Q  to  mgg^  can 
be  found 


log 


*"400 

"’254 


400 


'Jm  ‘'’8  m 


2.62(.  197)  .  .517. 


"’400 

*"254 


3.29. 


*"254  °  *"400* 


The  relationship  between  m  at  400  MHz  and  m  at  254  MHz  is  plotted  in  Figure  12. 
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We  have  already  related  the  AFCRL  SI  to  m  as  shown  In  Figure  7.  This  was 
based  on  experimental  '<ta  at  137  MHz.  Unless  the  scatter  mechanism  changes 
at  other  frequencies.  It  Is  reasonable  to  asstune  that  this  relationship  is  valid  at 
?.54  MHz.  Therefore,  from  Figures  5,  7,  and  12,  we  can  construct  the  desired 
relationship,  MHzl  vs  S^%(400  MHzi  by  the  following  steps: 

(a)  given  a  value  of  S^(400).  read  m  from  Figure  5; 

ex,  ®  10%,  hi^QQ  ®  100. 

(b)  given  read  from  Figure  12; 

ex.  •  100,  *"254  ’ 

(c)  given  m,  read  SI%  from  Figure  7; 

ex,  m  •  30,  SI  >  35%, 

Therefore,  an  3  10%  at  400  MHz  corresponds  to  an  SI  •  35%  at  254  MHz,  These 
steps  can  then  be  followed  to  plot  the  curve  shown  in  Figure  13. 
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9.  CONCLUSIONS 

Much  of  the  eaalyaie  that  hes  been  preaented  depends  on  a  knowledge  of  the 
frequency  dependence  or  spectral  Index  of  Ionospheric  scintillations.  A  constant 
power-law  function  of  frequency  was  assumed  for  scintillation  Index  and  m.  and 
in  the  100>  to  400-MIIb  range  this  may  be  an  error.  The  evidence  of  extreme 
scintillations  makes  this  assumption  particularly  tenuous  for  the  equatorial  regions. 
Additional  data  must  be  acquired  at  closely  spaced  frequencies  to  more  adequately 
define  the  dependence  of  scintillations  on  frequency. 
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Appendix  A 


httnietion  for  Completing  FORTRAN  Codm|  Form 
for  Scaling  15-m)n  Scintiilatlona  (ATS-3) 


Scintillations  are  scaled  from  the  chart  as  a  decibel  change  corresponding  to 

P  -  P  .  where  P  .  is  the  third  peak  down  from  the  maximum  and  P_,  is 
max  min  max  min 

the  third  null  up  from  the  lowest  excursion.  The  values  for  each  15*mln  period 
are  recorded  in  3  columns  (tens,  units,  and  tenths  of  dB>  starting  with  column  12. 
The  scintillations  are  measured  from  the  calibration  steps  as  a  relative  change 
in  decibels.  Each  card  is  a  4 -hour  block,  indexed  by  column  7,  Use  >99  to  record 
interference  or  inability  to  make  a  measurement  in  any  lS*mln  period.  Use  001 
to  record  signal  present  but  scintillations  too  sntall  to  read,  that  is,  scintillations 
less  than  1  dB.  When  there  is  interference  or  poor  data  during  an  entire  4 -hour 
block,  put  a  9  In  column  10.  When  signal  is  present  but  scintillations  are  too 
small  to  read  during  an  entire  4-hour  block,  put  a  1  in  column  lO;  the  remainder 
of  the  line  will  be  left  blank.  If  data  is  missing  for  a  full  4  -hour  period,  no  card 
is  required.  If  data  is  missing  or  poor  for  a  portion  of  a  4 -hour  period,  -99  is 
required  in  the  periinent  columns. 


Practiiif  Ml*  Unk 
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IS -MIN  SCINTILLATION  FORMAT 


1.2 

Year 

3.4 

Month 

s.e 

Day 

7 

4-hour  period  1  *  00  -  04 

2  •  04  -  08 

3  «  08  -  12 

4  •  12  -  18 

5  >  18  -  20 

8-20-24 

8.» 

Blank 

10 

9  >  Interference  for  4  hours 

1  ■  actntillatlona  less  than  1  dB  for  4  hours 

11 

Blank 

12, 13.  14 

lat  15 -min  scintillation  reading  in  dB 

12  *  tens  dB.  13  ■  units  dB.  14  ■  tenths  dB 

15,16. 17 

2nd  15 -min  reading 

etc.,  to 

etc.,  to 

57,58,  SO 

12th  IS -min  reading 

60-67 

Blank 

68,69 

Sources  A1  •  ATS-l,  AS  •  ATS-3,  A5  *  ATS 
12  *  Intelsat  2  P-2 

70,71,72 

Frequency 

73.74 

Station  No. 

70,80 

dB 

